^®A1 and the Formation of the Solar System from a Molecular 
Cloud Contaminated by Wolf-Rayet Winds 

Eric Gaidos^ 

Department of Geology and Geophysics, University of Hawaii, Honolulu, HI 96822 
' gaidos@hawaii.edu 

(N : 

d '. Alexander N. Krot^ 

' ' ' Hawaii Institute of Geophysics and Planetology, University of Hawaii, Honolulu, HI 96822 
T— I \ 

, sasha@higp.hawaii.edu 

^ . Jonathan P. Williams^ 

• , Institute for Astronomy, University of Hawaii, Honolulu, HI 96822 
^ ^, jpw@ifa.hawaii.edu 
Q . and 

, Sean N. Raymond^ 

NASA Postdoctoral Program Fellow, Center for Astrophysics and Space Astronomy, University of 
^ ■ Colorado, Boulder, CO 80309 

' raymond@lasp.colorado.edu 

\o : 

m ■ 
cn ; 

^ ■ 

O ; ABSTRACT 

. In agreement with previous work, we show that the presence of the short-hved radionuchde 

^ I ^^Al in the early Solar System was unlikely (< 2% a priori probability) to be the result of direct 

•'H , introduction of supernova ejecta into the gaseous disk during the Class II stage of protosolar 

' evolution. We also show that Bondi-Hoyle accretion of any contaminated residual gas from 

^ . the Sun's natal star cluster contributed negligible ^®A1 to the primordial Solar System. Our 

■ " " ' calculations are consistent with the absence of the oxygen isotopic signature expected with any 

late introduction of supernova ejecta into the protoplanetary disk. Instead, the presence of ^^Al 
in the oldest Solar System solids (calcium-aluminum-rich inclusions or CAIs) and its apparent 
uniform distribution with the inferred canonical ^^Al/^^Al ratio of (4.5 — 5) x 10^^ support 
the inheritance of ^^Al from the Sun's parent giant molecular cloud. We propose that this 
radionuclide originated in a prior generation of massive stars that formed in the same molecular 
cloud and contaminated that cloud by Wolf-Rayet winds. We calculated the Galactic distribution 
of ^^Al/^^Al ratios that arise from such contamination using the established embedded cluster 
mass and stellar initial mass functions, published nucleosynthetic yields from the winds of massive 
stars, and by assuming rapid and uniform mixing into the cloud. Although our model predicts 
that the majority of stellar systems contain no ^^Al from massive stars, and that the a priori 
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probability that the ^^Al/^'^Al ratio will reach or exceed the canonical Solar System value is only 
~6%, the maximum in the distribution of non-zero values is close to the canonical ^^Al/^^Al ratio. 
We find that the Sun most likely formed 4-5 million years (Myr) after the massive stars that were 
the source of ^^Al. Furthermore, our model can explain the initial Solar System abundance of 
a second, co-occuring short-lived radionuclide, ''^^Ca, if ~ 5 x 10^ yr elapsed between ejection of 
the radionuclides and the formation of CAIs. The presence of a third radionuclide, ^°Fe, can be 
quantitatively explained if (a) the Sun formed immediately after the first supernovae from the 
earlier generation of stars; (b) only 5% of supernova ejecta was incorporated into the molecular 
cloud , or (c) the radionuclide originated in an even earlier generation of stars whose contributions 
to other radionuclides with a shorter half-life had completely decayed. 

Subject headings: planetary systems: protoplanetary disks — planetary systems: formation 



1. Introduction 

1.1. Short-lived radionuclides in the early 
Solar System 

Primitive meteoritic materials contain com- 
pelling evidence for the short-lived radionucides 
(SLRs) lOfie, 26A1, 36ci^ 4iCa, ss^n, eopg^ 
^°^ Pd, and ^^^Hf in the early Solar System 



(jGoswami et al. 20051 ). These radionuclides have 
a half-life ti/2 < 10 Myr and are potential 
high-resolution chronometers o f events during 
the epoch of planet formation (jKita et al. 2005t 
Krot et al. 200 8a'). but only if they were in- 
troduced at discrete times and were unifor- 
maly distributed in the Solar System. The 
potential sources of these radionuclides inform 
us about the youn g Sun's stellar neighborhood 
(jHester et al. 20041) , or its magnetic i nteraction 
with an accretion disk (jShu et al. 19971 ). The de- 
cay of one short-lived radionuclide, ^^Al, might 
have been the principle heat source in plan- 
etesimals and responsible for the differentiation 
of the parent bodies of magmatic iron mete- 
orites in the first 1-2 M yr of the Solar System 
(iGreenwood et al. 2005t ISchersten et al. 2006t 
Markowski et al. 20061 ). 

The origin of SLRs is controversial. The half- 
life of each is much shorter than the '^ lOO Myr 



(WR) stars (jArnould fc Prantzos IQSi ): and (ii) 
production by irradiation of gas or dust with 
energetic particles from the active young Sun 
( Heymann fc Dziczkaniec 1976[ ). An alterna- 
tive scenario invoking an origin in a nearby 
i ntermediate-mass asyrn ptotic giant branch star 
(jWasserburg et al. 19941 ) is generally discounted 
because of the very low probabil ity of such an 
encounter ( Kastner fc Myers 19941) . Each of the 
two scho ols of thought h as developed elaborate 
models (iLee et al. 1991 iGounelle et al. 2006t 



Ouellette et al. 200?! ) but neither has produced 



a comprehensive explanation of the origin, abui> 



dance, and distribution of all the SLRs oswa mi et al. 20051 : 
Gounelle et al. 20061 : [Puprat fc Tatischeff 200?! ). 



The origins of three SLRs seem unambiguous: 
(i) ^°Be (ti ^ 1.5 My r) inferred from excess 



( McKeegan et al. 2000t ) could have been produced 
by energetic particle irradiation but not by stel- 
lar nucleosynthesis. Although magnetic trapping 
of Galactic cosmic rays in the protosolar molec- 
ular cloud has been advanced as an alternative 



mixing time of the interstellar medium (|de Avillez fc 
and the excess abundances of at least five ra- 
dionuclides (i°Be, 36C1, 26 ^^Ca, and ^°Fe) 
require one or more "local" sourc es in addition t o 
the average Galactic background ( Jacobsen 2005[ ). 
Two principle scenarios emerged soon after the 
first reports of fossil SLRs in meteorites: {i) 
an origin in one or more neighboring massive 
stars, either Type II superno va (SN) progen- 
itors ( Gameron fc Truran 19771 ) or Wolf-Rayet 



explanation (Desch et al. 20041 ). it is inconsistent 

with variations in the infer r ed initial -'^"Be/^Be 

ratio ('McKeeean et al. 2000*; 'Sugiura et al. 200ll : 

Marhas_et. al^ 20.02; MacPhcrson et al. 20031 ). (ii) 

Excess ^^S correlated with the ratio ^^Gl/^'^S and 

— — attribttteble to the decay of ^^Cl (ri/2 ~ 0.3 Myr) 
MacLow 20021, , . , , ,^ ; 
was iupurtcd m sodalite, an alteration phase that 

replaced anorth ite in GAIs and chondrules from 

GV chondrites (|Lin et al. 20051 : IHsu et al. 20061 ). 

The inferred ^^Gl/^^Gl ratio at the time of so- 
dalite formation is 5 x 10~^. If the sodalite formed 
late (>1.5 Myr after GAI crystallization, based 
on absence of ^SAl), the initial ^'^Cl/^^Cl ratio 
was > 1.6 X 10~*. This is inconsistent with a 
massive stellar source and requires a late episode 
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of irradiation. (Hi) ^°Ni excess correlated with 
the ^^Fe/^®Ni ratio is evidence for Uve ^°Fe, a 
radionuchde that cannot be produced by irradi- 
ation and must have originated in one or more 
massive stars ( Tachibana fc Huss 20031 ). High- 
precision nickel isotope measurements in several 
groups of magmatic iron meteorites indicate that 
^"Fe was uniformly distri buted in the solar neb- 
ula (jPauphas et al. 20081 ) but its initial abun- 



is uncer- 



danc e. r°Fer^Fe)n ~ (0.5 - 1) x 10 

tain dXachibana et al. 2006 : Bizzarro et al. 2007 



Quitte et al. 20071 : iGuan et al. 20071 ). The lower 
end of the range of estimates is consistent with 
the expected abundance in star-forming regions 
if star formation rates were approximately twice 
as hi gh at the epoch of Solar System forma - 
tion (jWilhams 2008t iGounelle fc Meibom 2008al ) 
and/or the half-life of ^°Fe is actually longer than 
the published value {ti/2 ~ 1.5 Myr). 

1.2. The origin of ^^Al 

Ironically, the origin of ^'^Al (t i ~ 0.73 Myr) 
the first SLR to be discovered (|Lee et al. 19761 ) 
and the one best studied, remains an enigma. A 
significant contribution by irradiation is disputed 
on the grounds that (i) production models ad- 
justed to achieve the required ^^Al/^^Al ratio over- 
predict the observed initial abundance of ''^Ca 

{ti/2 ~ 0.1 My r) whose co-occurrence ind i cates a 

common origin ( Sahiipal fc Goswami IQQSt Goswami et al 
(ii) the flux of energetic particle inferred from 
X-ray observations of young solar-type stars 
would have been insufficien t to produce the ^^Al 
(|Duprat fc Tatischeff 20071 ). and (m) a lack of cor- 
relation between the presence of ^^Al and ^°Be 



of ^^Al in massive stars, the mechanism and 
timing of its delivery to the early Solar Sys- 
tem remains unclear. Two models have been 
proposed: instability-induced mixing of gaseous 
SN ejecta during the mol ecular cloud core phase 
(jCameron fc Truran 19771 ). and injection of SLR- 
beari ng dust grains into th e later protoplanetary 
disk (jOuellette et al. 20051 ). The relatively small 
cross-section of the solar nebula dictated that the 
proge nitor was within ^ 1 pc of the Solar Sys- 
tem (jLoonev et al. 20061 ). Such a circumstance 
is possible only in t he dense environm ent of a 
large stellar cluster (jHester et al. 20041 ). Incor- 
poration of hot, low-density gas from SN ejecta 
into the denser, coo l er disk gas is ineff i cient 
(jVanhala fc Boss 20021: [Ouellette et al. 2007tlBoss 200S 



Instead, Ouellette et al. (2005) proposed that 
several SLRs, including ^^Al, were delivered 
to the early Solar System as grains that con- 
densed from SN ejecta and vaporized upon en- 
tering the relatively h igh-density gas in the disk 
(jOuellette et al. 20071) . Disks around low-mass 



cluster stars persist for up t o 6 Myr (jHaisch et al. 2005 
Javawardhana et al. 2006h . longer than the main 
seque nce lifetime of the m ost massive SN progen- 
itors (|Schaller et al. 19921 ). 

However, any explanation for Solar System 
^^Al invoking a late introduction of SN ejecta 
has four significant shortcomings : (?) Not 
a^QQ^ rs form in large clusters and most clus 



which was formed by irrad iation (jMacPherson et al. 2003 
Marhas fc Goswami 2004h . Finally, the canonical 
value [(4.5-5) x 10"^] of the initial ^^Al/^^Al ratio 
i n the majority of CAIs f ro m primitive chondrites 
(iMacPherson et al 19"9l IJacobsen et al. 2008t 
Makide et al. 200a) . the consistent chronology 
of CAI and chondrule formation between ^'^^Pb- 



ters are d ynamically unboun d and disperse in 
~10 Myr dlada fc Lada 20031 ). It is statistically 
unlikely that the Sun would have been suffiently 
close to a massive sta r at the end of the lat- 
ter's mai n sequence life (Williams fc Gaidos 20071 : 



^°^Pb and ^^Al-^^Mg sys te matics (lAmelin et al. 2002 ; 



Connelly et al. 20081 ). 



Hallidav fc Kleine 20061 : 
and the apparently uniform Mg-isotope compo- 
sitions of bulk chondrites, Mars, Moon and the 



Earth (jThrane et al. 20061 ) are evidence for a uni- 
form distribution thought inconsistent with a cen- 
tral (solar) irradiation source. 

Although these observations favor an origin 



Gounelle fc Meibom 2008bl ). (ii) Even the most 
massive stars have a main sequence lifetime of 
at least 3 Myr and, if they and the Sun formed 
simultaneously, the former would have ended in 
supernovae late in the evolution of the solar 
nebula and probably long after CAIs contain- 
ing the canonical ^®A1/^^A1 ratio had formed. 
These CAIs, the oldest dated so lids from the So- 
lar System ( Amelin et al. 20"o3 ). have a narrow 
range of inferred initial ^^Al/^ ^Al ratios suggesting 
that they formed in < Cl Myr (jThrane et al. 2006t 
Jacobsen et al. 20081 ). consistent wit h the dura- 
tion of Class O-I stages of pro tostars ([Smith 2004 : 
Ward-Thompson et al. 20071 ). A short interval of 
CAI formation is also consistent with the nar- 
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row range of their oxy gen isotope compositions 
(A^^O = -24 ± 2 %o) iMacPherson et al. 20081 : 



Palacios et al. 2005HSahiipal fc Soni 20061 ). 



Makide et al. 2008h . which are similar to the in- 



ferred oxygen isotopic compos ition of the Sun 
( Hashizume fc Chaussidon 20051) . [Later CAI 
formation would presumably have reflected the 
rapid oxygen isotopic evolution of the solar neb- 
ula along the slope-one carbonaceous chondrite 
anhydrous minerals (CCAM) line towards the 
terrestrial value (A^''0 = 0%c), a trend at- 
tributed to CO photochem ical self-shielding and 
radial mixing in the disk (jYurimoto et al. 2007 : 
Aleon et al. 2007t ).] (Hi) Late injection of ^""Al 



into the protoplanetary disk would likely have 
disturbed the oxygen isotopic composition of the 
disk from the CCAM line, a nd this is not observed 
(|Gounelle fc Meibom 20071 ). (iv) SN models that 
produce the canonical ^^Al/^^Al ratio invariably 
over-predict the abundance of ^'^Mn and must im- 
pose fallback of the innermost layers, e.g., Meyer 
(2005). (These models also over-predict the abun- 
dance of sope, see ^Oj) . 

A scenario in which SLRs from coeval SN pro- 
genitors were injected into the protoplanetary 
disk can be evaluated for its statistical plau- 
sibility. Gounelle fc Meibom (2007) estimated 
that the probability that any given disk is con- 
taminated by a SN with enough ^^Al to reach 
the canonical ^^Al/^'^Al ratio is < 0.3%. How- 
ever, there are three limitations to their model: 
(z) It underestimated the probability by assum- 
ing a maximum disk radius of 50 AU based 
on observations and theoretical arguments that 
disks within 0.1 pc of massive (O) stars suf- 
fer from photoevaporation (Johnsto ne et al. 19981 : 
[chevalier 200Q). But dis ks further from O 
stars are invariably larger ( Vicente fc Alves 2005t 
Andrews fc Williams 2007ilBalog et al. 2007[) and 



this effect compensates the greater distance from 
the source of radionuclides, {ii) Their model over- 
estimated the probability by neglecting the expan- 
sion of the host star cluster and using the <1 Myr- 
old Orion nebular cluster as a template, rather 
than a cluster at the minimum age (3 Myr) when 
supernovae occur. In older clusters, stars are more 
widely separated, (iii) It only included the contri- 
bution of SN to ^^Al. ^^Al is also produced and 
ejected from stars with initial masses > 40 Mq in 
their Luminous Blue Variable (LBV) and Wolf- 
Rayet (WR) phases ( Arnould fc Prantzos 1986t 



1.3. Was Solar System ^''Al inherited? 

The uniform distribution of ^^Al, the exis- 
tence of the canonical ^^Al/^^Al ratio in CAIs, 
the adherence of primitive oxygen isotope com- 
positions to the CCAM line, and the low like- 
lihood of a SN injection event all point to the 
introduction of ^^Al before the collapse of the 
protosolar cloud. Indeed, the oxygen isotopic 
composition of the entire Solar System appears 
to be disp laced from the loc us of mean Galactic 
evolution ( Young et al. 2008h . suggesting primor- 
dial contamination. The source of ^^Al must also 
have introduced ^^Ca, which co-occurs with an in- 
ferred initial abundance of ^ ^Ca/^°Ca ~ 1.5 x 10~^ 
(jSahiipal fc Goswami 19981 ). The half-life of this 
radionuclide is only 0.1 Myr, limiting any time 
delay between production at the source and its 
incorporation into CAIs. The source of ^^Al can- 
not have been accompanied by substantial ^^Mn, 
as is the case of SN eject a without fallback onto 
the remnant (jMever 20051 ). (The predicted ^■^Mn 
abundance in the interstellar medium is consistent 
with its inferred initial abundance in the Solar 
System and a "local" source is not required.) Any 
relationship between the source of ^^Al and that of 
^°Fe remains to be determined. There is as yet no 
evidence that ^°Fe and ^^Al are correlated and had 
the same origin. Indeed, there may be a deficit of 
^°Fe relative to ^^Al compared to SN ejecta that 
cannot be explained by free decay of the two ra- 
dionuclides: The initial ra tio of ^°Fe to ^^Al in th e 
Solar System was 0.1-0.2 (|Tachibana et al. 2006h . 
lower than the 0.3 deduced from the Galactic 



average 7-ray emission (iWang et al. 20071) and 
theoretical predictions ( Limongi fc Chieffi 20061 : 



Wooslev fc Heger 20071 ) . 

Very massive (> 40 Mq) 



stars eject ^^Al 



(and other SLRs) during the Wolf-Rayet phase 
of mass loss near the end of hydrogen core burn - 
ing, as well as in SN ( Arnould fc Prantzos 19861 ). 
WR winds might account for a large fraction of 
the total fluence and Galac tic distribution of 7- 
rays from the decay of ^^Al ( Palacios et al. 20051 : 
iDiehl 2006t IVoss et al. 20081: iMartin et al. 2008^ . 
[The non-detection of 7-ray emission in the de- 
cay line of ^^Al from the nearest WR star 7^- 
Velorum can be explained by the dispersal of 
most of the radionuclide to large angular separa- 
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tion (|Mowlavi fc Mevnet 20061 )]. We propose that 
most or all of the ^^Al in the early Solar System 
originated in WR winds from one or more massive 
stars that contaminated the molecular cloud from 
which the Sun formed. (Sahijpal & Soni (2006) 
also considered the contribution of WR winds to 
Solar System inventories of SLRs.) These stars 
could have been members of the same embed- 
ded cluster as the Sun, or, more likely, members 
of another cluster that formed in the same giant 
molecular cloud (GMC) (Figure [1]). An analogous 
"self-contamination" scenario has been invoked to 
explain anomalo us abundanc e patterns in some 
globular clusters ( Smith 2006f l. 

Our proposal is based upon the following: (i) 
The amount of Al ejected in the win ds of a single 
60 Mq star (jLimongi fc ChiefH 200^ is sufficient 
to contaminate 2 x 10"* Mq of solar-metallicity 
gas to the canonical ^^Al/^^Al ratio of the Solar 
System, (ii) The WR phase can occur as soon as 
1-2 Myr after st ar formation and mor e than 1 Myr 
before the SN (jPalacios et al. 20051 ). making it 
more likely to contaminate residual star-forming 
molecular gas. (i n) WR winds have speeds of up to 
2000 km sec-i (jNiedzielski fc Skorzvnnski 20021 ) 
and can traverse star- forming regions (10-100 pc) 
in 10^ — 10^ yr. (iv) In contrast to single clusters 
where star-formation may be co-eval, star forma- 
tion in a molecular cloud can occur over an inter- 
val of at least a few My r (Hartmann et al. 2001), 
and possibly longer ( Williams et al. 2000l ) as 
clumps o f gas with a mass spectrum M~-^-'^^^-^ 
(jPudritz 2002 ) form multiple embedded stellar 
clusters (fwilliams et al. 2000l ). For example, 
the Orion star-forming complex contains sev- 
eral subgroups that are se veral Myr o lder than 
the Orion nebula cluster (|Ballv 200^ (v) Wolf- 
Rayet winds contain multiple SLRs, including 
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Al, -^^Ca, and ^^Cl, but httle or no ^^Mn or 
6°Fe (jArnould et al. 2006h . 



We propose that 
the collapse of the protosolar cloud homogenized 
the distribution of these isotopes (but see H5.2p . 
Our scen ario does not preclude a SN-triggere d 
collapse (|Cameron fc Truran 19771 IBoss 2008f) . 
which would have occured after the WR phase. 

In f2]we revisit the scenario of Al-introduction 
by SN into the protoplanetary disk with a Monte 
Carlo approach that used more realistic disk sizes, 
accounted for the expansion of clusters, and in- 
cluded the contribution from both WR winds and 



SN. While our results differ quantitatively from 
those of Gounelle fc Meibom (2007), the calcu- 
lated probability of a disk having the canonical 
Solar System ^^Al/^^Al ratio is nevertheless small 
(<2%). We also consider Bondi-Hoyle accretion 
of contaminated intracluster gas, as proposed by 
Throop fc Bally (2008) and find its inclusion does 
not significantly alter this result (Q. We then use 
a similar Monte Carlo model to investigate a sce- 
nario where ^^Al is introduced by WR winds into 
the parent molecular cloud of the Sun (Q. This 
model readily reproduces the canonical ^®A1/^^A1 
ratio of the Solar System. We discuss the delivery 
of ^^Al from WR winds into the molecular cloud, 
interpret CAIs that lack the canonical ^^A1/^''A1 
ratio, and, in the context of the cloud contamina- 
tion scenario, present additional calculations for 
three other SLRs (^^Ca, ^Ope, and ^^Cl) in gl 

2. Disk contamination scenario 

We calculated the Galactic distribution of 
^^Al/^^Al ratio in t he disk injection scenario 
(|Ouellette et al. 2005h by extending the model 
of WiUiams & Gaidos (2007) to include '^^Al from 
both winds and SN, a dynamically realistic de- 
scription of cluster expansion, and the effect of 
the UV field of massive stars on disk size. Un- 
like Williams & Gaidos (2007) we do not explic- 
itly consider disk evolution and disappearance in 
our mode l because this occu rs on a timescale of 
~6 Myr (|Haisch et al. 2005f ). much longer than 
the likely epoch of CAT formation. 

For simplicity, we assume that every young 
low-mass star has a disk and that the properties 
of the disk are independent of stellar mass. Ac- 
cording to Lada & Lada (2003), 72% of stars form 
in clusters with N,, > 100 members distributed 
in size according to dNc/dN^, oc N^"^. The re- 
maining 28% form in isolation or in clusters with 
fewer than TV* = 100 stars. The most massive 
members of such small clusters will be ~5 Mq B 
stars that do not produce winds or Type II core- 
collapse SN. Disks around members of such clus- 
ters will not receive any exogenous ^^Al. The 
other 72% of disks were represented by 10^ Monte 
Carlo calculations. The size of the host cluster 
of each disk was drawn from a A^^^ distribution 
(A* > 100). The number of massive stars (SN pro- 
genitors with Af* > 8M0) in the host cluster was 
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selected from a Po isson distribution with an aver- 
age of 3 X 10-3 AT* (I Williams fc Gaidos 20n7h . The 



masses of these stars were drawn from a power- 
law initial m ass function dN^,/dM^, (x M~'^-^ 
( Kroupa 2002 ) . Clust ers form over ~1 Myr 
([Hille nbran d et al. 20071 ) and most star forma- 
tion in a single cluster occurs within ^3 Myr 



(iHartmann 200lVHartmann 2 003trHuff fc Stabler 20061 : 
[Jeffries 2007: Hillenbrand ct al. 2007[ ). We as- 
sumed the instantaneous formation of all massive 
stars and an exponentially decaying rate of low- 
mass star formation after massive star formation. 
We used age sta tistics for members o f the Orion 
Nebula Cluster ([Palla fc Stabler 1999[ ) to infer a 
decay time of 2.7 Myr, which we adopted for clus- 
ters of all sizes. We assumed that the rate of star 
formation does not depend on stellar mass. 

The ^^Al/^'^'Al ratio in the disk around the i- 
th low-mass star at an interval T after the star's 
formation at time t* was calculated by summing 
over the product of the the yield rrij of the ra- 
dionuclide from the j-th wind or SN, the solid an- 
gle subtended by the disk at the time of ejection 
tj , and the factor of free decay between that time 
and t* + T- 

(1) 

where Vd and rrid are the (constant) radius and 
mass of the disk, Jm is the mass fraction of ^^Al, 
and dij is the distance of the disk from the source 
star during the at the time of ejection. We con- 
sidered injection as instantaneous because the 
speeds of SN ejecta and winds are > 1000 km s"^, 
and the ^^Al-producing WR phase of an individ- 
ual massive star is brief (^^1 Myr) compared to 
the dispersal timescale of the cluster (~10 Myr). 
To account for the changing perspective of each 
star as it orbits inside the cluster, we used the 
isotropic average of the projected cross-section 
of each disk, 7rr^/2. rrij and tj were estimated 
by spline interpolation in a grid of yield calcula- 
tions by Limongi fc Chieffi (2006). We used a de- 
fault disk radius of 200 AU ([Vicente fc Alves 2005t 



Andrews fc Williams 20071 ) but to account for 



photoevap oration by the UV rad iation from mas- 
sive stars ( Johnstone et al. 19981 ) we reduced this 
to 30 AU for disks within 0.2 pc of the cluster 
center at 3 Myr, and to zero for disks within 



0.1 pc. This was a conservative assumption, since 
the evidence for sig nificant disk truncation is weak 
([Balog et al. 20071 ). but it only has a minor effect 
on our results. Like Ouellette et al. (2005), we 
adopted a Minimum Mass Solar Nebula disk mass 
of 0.013 M0 and the fractional abundance of ^^Al 
given by Lodders (2003). 

As in Williams fc Gaidos (2007), we placed the 
massive stars at the cluster center ([Crebel 20071 ) 
and assumed that, at any time, low-mass stars 
were distributed with an inverse-square density 
profile, such that an equal number of stars reside 
in shells of constant thickness out to the edge of 
the cluster at Tc- To model cluster expansion, and 
thereby determine the distance of a disk to each 
source of Al, we developed an empirical relation- 
ship for the time dependence of Vc based on a se- 
ries of numerical simulations of clusters containing 
between 100 and 15,000 stars. The dynamical sim- 
ulations were preformed using the NB0DY4 code 
runni ng on the Cam bridge University GRAPE-6a 
card ( Aarseth 20031 ). In each case the stars were 
initially (3 Myr) distribute d in a Plummer sphere 
([Binnev fc Tremaine 19881 ) with a radius set by 
the requirement that the initial surface density 
S3 = 100 pc-2 ([Williams fc Gaidos 20071) . The 
virial parameter VL (ratio of kinetic to gravitational 
potential energy) was set to 1.5. This condition is 
brought about by a cluster local star formation ef- 
ficiency of 33% in the embedded cluster and the in- 
stantaneous removal of the remaining gas at 3 Myr 
( Bastian fc Goodwin 20061 ). The size of the clus- 
ter was explicitly determined at regular intervals 
until an age of 10 Myr. We found that the ex- 
pansion of the cluster from its size at 3 Myr was 
closely approximated by 



rc(t) «r,(3Myr)+0.45 



2GM*f7V7i^VV* 



1/2 



(i-3Myr), 
(2) 

where G is the gravitational constant and Af* 
is the average stellar mass. The expression 
multiplying 0.45 is the cluster's virial speed 
([Binnev fc Tremaine 19881 ). In each of our sim- 
ulations the surface density of the cluster fell to 
the background level of field stars (2-3 pc"^) 
by 10 Myr, in a greement with observations 
([Lada fc Lada 20031 ). 

Eqn. [2] specifies the radius of the cluster at 
the epoch tj of a massive stellar wind or SN. The 
uniform distribution of low-mass stars with dis- 
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tance from the cluster center < dij < rdtj) was 
then used in Eqn. [T] to produce a distribution of 
^^Al/^^Al. These distributions were summed over 
all events in a cluster, corrected for free decay, av- 
eraged over 10^ realizations, and multipled by 0.78 
to produce a Galactic distribution. 

The calculated distributions of ^^Al/^^Al are 
plotted in Figure [5] and the fractions of systems 
that have any initial ^^Al or ^^Al abundances 
exceeding the canonical value are given in Ta- 
ble [TJ We also report the 95-percentile values 
of ^^Al/^'^'Al. We considered two values for T, 
which in the Solar System represents the epoch 
of CAI formation. The probability is 1.1% for 
T = 0.5 Myr, rising to 1.9% by 1 Myr. The 
probability is higher at still later, but unlikely 
CAI formation times (not shown). Our probabli- 
ties are several times higher than that reported 
by Gounelle & Meibom (2008), mostly due to the 
larger disk size we used, but are still small. Our 
estimates are nonetheless optimistic because we 
assume all high-mass stars form prior to low-mass 
stars. If all stars form instantaneously, then nearly 
3 Myr must elapse before ^^Al is produced, and no 
disks contain ^^Al by the time of CAI formation. 

3. Bondi-Hoyle accretion scenario 

Disks might also accrete gas from their na- 
tal cluster before it is removed by winds and 
SN explosions. Throop & Bally (2008) pro- 
posed that Bondi-Hoyle accretion of residual, 
SN-contaminated cluster gas onto the protosolar 
nebula produced isotopic anomalies in the Solar 
System, including the presence of SLRs. They 
estimated that disks around solar-mass stars in 
a 3000-star cluster could accrete an additional 
^0.01 M0 of gas, an amount comparable to the 
original mass of a disk, over the 2-4 Myr that 
gas remained in the cluster. We estimated the 
amount of ^^Al that could be introduced by this 
process by specifying the fraction of disk mass / 
acquired by Bondi-Hoyle accretion by the time 
of CAI formation, and assuming that accretion is 
constant as long as intracluster gas is present. We 
assumed that the maximum amount of mass that 
can be acc reted by a disk is equ al to the initial 
disk mass (iThroop fc Ballv 2008f) . i.e. / < 0.5. 
To account for the absence of cluster gas during 
the history of later-forming stars we adjusted the 



accreted mass by the ratio of the time interval be- 
tween low-mass star formation and the disappear- 
ance of cluster gas, and the lifetime of the cluster 
gas. Like Throop & Bally (2008), we assumed a 
gas lifetime of 3 Myr. We calculated the average 
^^A1/^''A1 ratio of cluster gas during the period 
of Bondi-Hoyle accretion (see []4] for details) and 
then determined the final ^^A1/^''A1 of the disk as 

26A1/27A1 = (l-f)(26Al/27Al)o + f(26Al/27Al),ec 

where the subscripts refer to the initial value and 
the average value during accretion. 

We calculated distributions of the ^^A1/^''A1 ra- 
tio for combined Bondi-Hoyle accretion and disk 
injection and found that the former has a negligi- 
ble effect (Figure [2] and Table [T]). This is because 
low-mass stars that form late enough (~3 Myr) to 
acquire significant ^^Al will accrete little gas be- 
cause the intracluster gas disappears soon there- 
after. This is of course entirely a result of our (rea- 
sonable) assumption that intracluster gas is evac- 
uated by the time the massive stars leave the main 
sequence, if not earlier. Nevertheless the same cos- 
mochemical timing arguments that apply to the 
disk injection scenario also apply to the Bondi- 
Hoyle scenario; CAIs probably formed by the time 
a disk had formed, or very soon thereafter, and 
thus later accretion of contaminated gas cannot 
be responsible for the presence of ^^Al in them. 

4. Molecular cloud contamination sce- 
nario 

In this scenario the Sun's natal giant molecu- 
lar cloud spawned an earlier generation of mas- 
sive stars [Figure [TJa)] whose WR winds contami- 
nated the rest of the cloud (b) , from which the Sun 
subsequently formed, perhaps in a second clus- 
ter (c and d). We calculated the distribution of 
26^j^27^j ratios in this scenario by 10*^ Monte 
Carlo simulations, each corresponding to a disk 
formed from a CMC contaminated by an immedi- 
ately previous generation of massive cluster stars. 
The number and masses of those stars were drawn 
from the distributions described in Sj2j The mass 
of gas in the CMC was calculated using an average 
stellar mass derived from the initial mass function 
of Kroupa (20 02), and a total star- formation effi- 
ciency of 10% (| Williams et al. 2000l ). The amount 
of ^^Al added to the molecular cloud by massive 
stellar winds was calculated as a function of time 
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using the yields and times of Limongi & ChiefB 
(2006) and by assuming 100% delivery efficiency 
and instantaneous mixing into the cloud (we dis- 
cuss this assumption further in ^5.ip . We did not 
include SN ejecta in these calculations, but include 
it when we consider ^°Fe in i|5.3l (but see the foot- 
note on SN ejecta delivery in §5.ip . We assumed 
that the Sun and CAIs formed simultaneously. 

We first carried out a series of calculations for 
different intervals of elapsed time (3-6 Myr) be- 
tween the formation of the earlier generation of 
massive stars and the Sun (dashed lines in Fig- 
ure [3]). If the interval is less than 3 Myr contam- 
ination has yet to take place in our model and 
newly-formed stellar systems lack ^^Al. An in- 
terval of 4-5 Myr is most likely to produce the 
CAI value. The history of low-mass star forma- 
tion in molecular cloud co mplexes is poorly known 



but clearly non-m onotonic (jHartmann et al. 2001 



Hartmann 2003h . We calculated more realistic dis 



tributions by again adopting the exponential rate 
with a 2.7 Myr decay time based on the data of 
Palla & Stabler (1999). This is plotted as the 
heavy solid line in Figure [3l That curve can be 
understood as a convolution of the star formation 
history with the ^^Al/^^Al distributions for "star- 
burst" scenarios. The probabihty of ^^A1/^''A1 
exceeding the canonical value is 6.2%. This fig- 
ure depends on the assumed star formation his- 
tory: For example, varying the decay time con- 
stant by ±1 Myr changes the fraction between 4.3 
and 6.5%. However, the peaks in all three distribu- 
tions are near the CAI value (solid lines in Figure 
[3]). This robustness is a result of negligible ^^Al 
production at times earlier than 3 Myr and negli- 
gible low-mass star formation at times much later 
than 5 Myr. We found that simulations which pro- 
duce a ^^Al/^^Al ratio within 2cr [a = 0.1 x 10"^ 
(jGoswami et al. 20051 )] of the canonical value were 
most likely to involve contamination of a GMC 
having ~ 3 x 10^ Mq of gas by the massive mem- 
bers of a cluster with ~ 10^ stars. Such a situ- 
tation is exemplified in our galaxy by NGC 3603, 
which contains multi ple Wolf-R ayet stars within 
an HII region (Melena et al. 200a) . 



5. Discussion 

5.1. Delivery of ^^Al into the molecular 
cloud 

Our scenario requires a plausible mechanism 
for the efficient introduction of the ^^Al carrier 
in Wolf-Rayet winds into the surrounding molec- 
ular cloud. In the calculations above, we as- 
sumed a 100% delivery efficiency but this will 
clearly not be the case. In general, mixing be- 
tween the hot, tenuous gas from massive stars 
and cooler, denser molecular gas is thought t o 
be very inefficient ( de Avillez fc Mac Low 2002[ ). 
SN-enriched gas from HII regions is thought to 
find its way back into the interstellar medium (if 
at all) through a circ uitous route taking '^100 Myr 
(|Tenorio-Tagle 2000l ). . mpact of SN ejecta onto 
a protostellar cloud core may induce mixing via 
Rayl eigh- Taylor and Ke l vin-Helmho ltz instabili- 
ties (|Boss fc Foster 1998t IBoss 200 8') but the in- 
jection efficiency is too low to explain the Solar 
System's canonical ^^A1/^''A1 raticQ. Gas in Wolf- 
Rayet winds will be less dense than SN ejecta by 
several orders of magnitude and efficient mixing in 
the gas phase is even less likely. Instead, the high- 
velocity (500-2000 km sec~^) winds will develop a 
reverse shock up on encountering the much denser 
molecular cloud (jWeaver et al. 19771 ). 

We propose that refractory dust grains are 
the principle carrier of ^®A1 and ^^Ca in Wolf- 
Rayet winds and that these can dynamically de- 
couple from the shocked wind and imbed them- 
selves into the surrounding molecular cloud, anal- 
ogous to the "aerogel" model described by Ouel- 
lette et al. (2005). Pre-solar grains of AI2O3 
(corundum or other forms) are found in primitive 
meteorites but their oxygen isotopes indicate a 
source in red giant or AGB stars, n ot in very 



massive stars (jHutcheon et al. 19941 lOtt 2007 : 
Nittler et al. 2008f) . Although many of these 
grains have large ^^Mg excesses produced by the 
decay of of ^^Al, the abundance of these grains is 



^ Adopting; an ^''Al mass fraction of 5.8 X 10 ^ 
l lLodders 20031 ). the protosolar cloud core initially 
contained 3 X 10"^ Mq of '^^Al. The mass fraction of ^^Al 
in the convective hydrogen shell of a 25 M p star at th e 
end of its main-sequence life is ~ 1 X 10-*^ iMcver 2005]). 
Therefore at least 3 X 10~^ Mq of SN ejecta must have 
been introduced (assuming no free decay). Boss et al. 
(2008) report that only 5 X 10~^ Mq is injected in their 
model. 
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insuffi cient to account for th e canonical ^^Al/^^Al 
ratio (jHutcheon et al. 1993 ). In fact, the ^""Al in 
CAIs must have been processed through the gas 
phase and subsequently recondensed into the re- 
fractory inclusions. The pre-solar grains are also 
much larger ('^l /im) than the silicate grains pre- 
dicted to form from winds and ejecta. Their size 
may be why they survived the incorporation pro- 
cess and the latter did not. 

Simple models of grain nuclcation and growth 
predict oxide grai n growth to sizes o f 0.01-0.1 
lim. in SN ejecta (jNozawa et al. 20071 ). SN are 
predicted to be copious sources of dust, but ob- 
servations have so far produced evidence only 
for a few tim es 10~^ Mp) of dust in individ- 



ual events (jMeikle et al. 20071) . Dust produc- 
tion in WR stars is poorly investigated, al- 
though such stars appear to be minor contrib- 
utors to the over- all interstellar dust budget 
(jTielens et al. 20051 ). Copious amorphous car- 
bon dust is observed around carbon-rich Wolf- 
Rayet CO stars and is thought to be the re- 
sult of colliding stellar winds in binary systems 
(jCrowther 2003[ ). WCO stars do noi produce ^SAl, 
but dust production by predecessor LBV and WN 
phases predicted to contain ^^Al in their winds has 
recently been establ ished (jRaiagopal et al. 2007 : 
iBarniske et al. 2008f ). The ry Carinae LBV star 
ejected ~10 Mf^^of dust-r ich materal during its 
1843 eruption (ISmith et al. 2003). including alu- 
minum oxide ( de Koter et al. 20051 ). although its 
^^Al content has yet to be de finitely establis hed 
by 7-ray observations (Knodlse der et al. 19961 ). 

To reach the molecular cloud, grains must sur- 
vive sputtering after they pass the reverse shock 
and move at high-speed 10'^ km sec~^) with re- 
spect to the shocked gas in the HII region. They 
also must not be completely decelerated within 
the reverse shock zone. These conditions place 
a lower limit on grain size. Grains too small 
are sputtered to destruction or eventually v apor- 
ized in the hot, shocked gas (jNozawa et al. 2007. ). 
The density of the wind 1 pc from a Wolf-Rayet 
star losing mass at a rate of 10~^ Mq yr~^ is 
10"^ cm~^ and the stopping distance of grains 
even as small as 0.01 fxm. grains is 10 pc, compara- 
ble or larger to the size of HII regions. The deceler- 
ation across the scale of the shocked region (~1 pc) 
(jWeaver et al. 19771 ) will be low and the fraction 
of material sputtered from the grains, which is re- 



lated to the deceleration [Eqn. 1 in Nozawa et al. 
(2007)] will be likewise smalB. 

Grains that escape the WR wind will not pen- 
etrate far into a GMC. Typical hydrogen num- 
ber densities in clouds are 10^ 



10'^ cm ^ 



and 



the stopping distance of 0.01-0.1 fim grains will 
be only of order ~ lO'^ AU. Gas densities in por- 
tions of the molecular cloud that are shocked and 
swept up by the expanding wind or SN ejecta will 
be higher and the stopping distances proportion- 
ally shorter. Thus, only the surfaces of clouds will 
be initially contaminated by ^^Al. Further trans- 
port of dust grains into cloud depends on their 
kinematics, which arc poorly understood. The 
smallest scale on which turbulence in clouds can 
affect the mass distribution and cause mixing is 
the sonic transition where the turbulent velocity 
is equal to the sound sp eed; this is rou ghly 1 pc in 
solar-metallicity clouds (jPadoan 1995f ) . Thus mix- 
ing might be inefficient at the surfaces of clouds. 
The degree to which this controls the incorpora- 
tion of ^^Al into new low-mass stars depends on 
the extent of large-scale mixing and whether star 
formation is triggered or at least assisted, by the 
i nteraction of winds or SN ejecta with cloud gas 
(jZavagno et al. 20071) . In that case, star forma- 
tion is spatially correlated with ^^Al abundance. 

Wolf-Rayet progenitors may themselves mi- 
grate into and contaminate regions of a molecular 
cloud where low-mass star formation has yet to 
occur. 10-30% of O stars have large peculiar veloc- 
ities (up to 200 km sec~^) relative to most early- 
type stars as a result either of dissolution of binary 
systems by SN expl osions, or encounters be tween 
two binary systems (jHoogerwerf et al. 200lh . One 
of these so-called "runaway" O stars moving at a 
typical speed of 30 km sec~^ would traverse a 
molecular cloud in ^1 Myr. Mass loss from this 
star, if occuring, could contaminate a larger region 
of the molecular cloud with SLRs. 



^In contrast, the density of shocked SN ejecta is ~ 10^ cm~^, 
the stopping distance of a 0.1 fim grain is only ~2 AU, 
whereas the ejecta scale length can be as large as 1 pc. 
Thus grains in SN ejecta are more likely to be trapped in 
the ejecta and never introduced into star-forming molecular 
gas. This is another argument for WR winds as the source 
of 2^A1 in the Solar System. 
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5.2. CAIs with low initial ^^A\p^A\ ratios 

Any scenario that explains the canonical Solar 
System ^^Al/^^Al ratio must also accommodate 
the exceptions. Several classes of CAIs show either 
no excess of ^^Mg produced by the decay of Al or 
have an inferred ^^A1/^''A1 ratio ^ 1 x 10^'^, much 
lower than the canonical value of (4.5-5) x 10^^. 
These include (i) igneous CAIs associated with 
chondrule-like materials (relict CAIs inside chon- 
drules and CAIs surrounded by chondrule-lik e, 



ferromagnesia n silicate rims) (jKrot et al. 2005a : 
Krot et al. 2005bl : {ii) some igneous CAIs in 
metal-rich (CB and C H ) carbonaceous cho ndrites 



([Gounelle et al. 2007t Krot et al. 2008b[ ): (m) 



FUN (fract ionation a nd unidentified nuclear ef- 
fects) CAIs (jLee 19881). iiv) isolated platy hi bonite 
crystals fPLACs) ( |lreland fc Feglev 2000l). iv] 



pyroxene- hibonite spherules ([Ireland et al. 1991 
Russell et al. 19981): ml some corundum-bearing 
CAIs (jSimon et al. 20021 ) . and {vii) most of the 



grossite- and hibonite-rich in cl usions in CH chon- 
drites (iKimura et al. 1995 IWeber et al. 1995t 
Krot et al. 2008bl ). 

{i and it): CAIs associated with chondrule 
materials and some igneous CAIs in CB and 
CH carbonaceous chondrites (Krot et al. 2001, 
2005a, b, 2008b) are 0l6-depleted to varying de- 
grees (A^^O ranges from -25%o to -5%c) relative 
to typical CAIs in primitive chondrites which 
unifor mly have ^^0-rich compositions (A^ '^0 ~ 
-25%o) (jltohet al. 2004 iMakide et al. 20081 ). We 
infer that the ^^Al-poor and 0l6-depleted CAIs 
experienced late-stage melting and oxygen isotope 
exchange, probably during chondrule formation, 
which could have reset their ^^Al-^^Mg systemat- 
ics. 

{iii-vii): The lower than the canonical ^^Al/^^Al 
ratio in FUN CAIs, PLACs, pyroxene-hibonite 
spherules, some corundum-bearing CAIs, and 
most of the grossite- and hibonite-rich inclusions 
in CH chondrites can be explained by (a) their 
late formation, after decay of ^^Al, (b) their early 
formation, prior to introduction of ^^Al, or (c) the 
lack of the canonical budget of ^^Al in their pre- 
cursors. Most of these CAIs h a ve 0l6-rich compo- 
sitions (ICoswami et al. 2001 : Simon et al. 20021 



canonical ^^Al/^^Al ratio. The rapid evolution of 
the oxygen isotopic composition of the inner So - 
lar System (|Krot et al. 2005d : lAleon et al. 20071 ) 
and the short (< 10^ yr) duration of CAI for - 
mation ( Thrane et al. 2006HJacobsen et al. 20081 ) 
thus make (a) unlikely. Similar arguments can be 
used against (b). Although an early formation 
is possible if complete melting and exchange of 
oxygen isotopes occured, this seems unlikely con- 
sidering evidence for inco mplete melting of CAIs 
(jMacPherson et al. 20051 ). 

We infer that the CAIs of categories iii-vii 
formed contemporaneously with ^^Al-rich inclu- 
sions. The absence of canonical ^^Al in their 
precursors suggests either they formed prior to 
homogenization of ^^Al in the Solar System 
teahiipal fc Gos wami 1995 iKrot et al. 2008bl: 
J<!rot et a l. 2008q ) or preferential loss of the (uni- 
formly distributed) ^^Al carrier d uring therma l 
processing of the CAI precursors ( Wood 1996[ ). 
Both explanations can be reconciled with the 
presence of nucleosynthetic a nomalies in some of 
these CAIs (|Lee et al. 19981) if the ^^Al carrier 
contributed a distinct component to the Solar Sys - 
tem's stable isotope composition (jLee et al. 19981 ). 
The second explanation is more speculative, how- 
ever. It hypothesises that (a) the precursors of 
these CAIs was isotopically heterogeneous and, 
contrary to typical refractory inclusions, escaped 
a cycle of complete evaporation-condensation, (b) 
the carrier of ^®A1 was relatively volatile, and (c) 
it was lost to varying degrees by sublimation of 
these CAI precursors prior to their melting. Al- 
though these inclusions can be used as an evidence 
for heterogeneous distribution of ^^Al among CAI 
precursors, the scale of any heterogeneity was 
probably limited because such inclusions are rare 
relative to ^""Al-rich CAIs. 

5.3. Other SLRs 

''-'^Ca: A further test of the wind model is 
whether it can also reproduce the inferred initial 



Krot et al. 2008bl : iKrot et al. 2008d) . indistin- 
guishable from those of typical CAIs with the 



■iiCa/^oCa ratio of 1.5x10"^ ()Sahiipal fc Goswami 199S 
Published calculations of ^^Ca yields in winds 
from massive stars are limited. We considered a 
60 Mq progenitor for which *^Ca and ^^Al yields 
from t he Wolf-Rayet winds were separately pub- 
lished (|Arnould et al. 20061 : iLi^ongi fc Chieffi 20061 ). 
We accounted for the additional free decay of Al, 
which is ejected during the hydrogen-burning WN 
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Wolf-Rayet phase, while ^^Ca is produced in the 
later core He-burning WCO phase 1.9 x 10^ yr 
later. The corrected ratio of "^^Ca to ^®A1 in the 
ejecta is ~ 25 times higher than in the Solar Sys- 
tem, but would be consistent if an additional time 
A « 0.5 Myr elapsed before CAI formation. It 
is interesting that this is approximately the same 
duration as the Wolf-Rayet phase itself before the 
final SN Ib/c. 

^''Fe: Wolf- Rayet winds contain negligible 
amounts of ^"^Fe (jArnould et al. 20061) . Live ™Fe 
in the early Solar System could have originated in 
SN from the same early generation of massive stars 
that produced t he ^^Al, or in an even earlie r gen- 
eration of stars (jGounelle fc Meibom 2008ai) . We 
repeated the calculations described in §4] but in- 
cluded SN contributions and calculated ^°Fe/^^Fe 
ratios in the same manner as ^^Al/^^Al, using the 
yields of Limongi & Chiefli (2006) and the solar 
iron abundance of Lodders (2003). We assume 
A — 0.5 Myr based on the ^^Ca abundance. In 
Figure m we plot Monte Carlo realizations for dif- 
ferent epochs (3-8 Myr before the Sun) for the ear- 
lier generation of massive stars. If all SN ejecta 
is incorporated, the ^'^Fe/^^Fe is over-predicted 
by an order of magnitude relative to ^^Al/^^Al. 
A comparision with the ratio of the two SLRs 
(black line) inferr ed from 7-ray measurements 
(jWang et al. 20071 ) suggests that this discrepancy 
may be in part the result of an overprediction of 
^°Fe yield - or underprediction of ^^Al yield - by 
the nucleosynthesis models (see Sj6]) . There are two 
other explanations suggested by Figure H) (a) the 
Sun formed 3 Myr after the massive stars, when 
many massive stars were in the WR phase but 
few SN had occurred (a scenario represented by 
the red dots extending below the primary locus); 
or (b) the SN contribution was attenuated by an 
effect such as described in §5.11 

Explanation (a) is statistically unlikely if star 
formation is uncorrelated and demands precise 
timing between the formation of massive stars and 
the Sun (^^3 Myr later), but could be demanded in 
a scenario where the Sun's formation was triggered 
by a SN (jCameron fc Truran 19771 : [Boss 2008) . 
Our simulations indicate that the initial solar 
26^j^27^j and ^°Fe/^^Fe ratios can be reproduced 
in this manner only in star clusters with > 10^ 
whose most massive members have ~ 100 M0. As 
such large clusters are relatively rare, this scenario 



is a priori less likely. Explanation (b) requires 
a reduction in the SN contribution to ^"Fe (and 
^^Al) by a factor of 20. This could be due to a 
combination of effects; retention or fallback o f the 
central region of the progenitor ( Mever 20051 ). in- 
efficient delivery of SN ejecta into the cloud ( ^5.1\i 
or the collapse of the protosolar cloud and a de- 
crease in its cross-section by the time the SN ejecta 
arrived. An alternative scenario (c) is that ^°Fe 
in the Solar System is the relict of an even ear- 
lier episode s of massive star-formation and con- 
tamination ( Gounelle fc Meibom 2008ah of which 
the ''^Ca and most ^^Al has decayed. ^°Fe will 
decay to 5% of its initial abundance in 6.5 Myr, 
during which ^^Al decays to 0.2%, and ''^Ca essen- 
tially vanishes. This last explanation is viable only 
if this earliest generation of stars ceased to con- 
tribute SLRs to the molecular cloud after ~6 Myr. 

•^^Cl: We estimated the abundance of ■^^Cl, 
which is also ejected by WR stars during the WCO 
phase (:Arnould et al. 2006[ ). and calculated the 
■^^Cl.'^^Cl ratio in the same manner as "'^Ca/'^'^Ca. 
We find that our model underpredicts the ratio by 
at least three orders of magnitude. It has already 
been recognized that stellar nucleosynthetic mod- 
els cannot account for this isotope, especially if it 
was introduced at the epoch of CAI formation 1- 
2 Myr before the host sodalite alteration phases 
were formed. At the present time, the only viable 
explanation appears to be a late episode of ir radi- 
ation by energetic particles ( Hsu et al. 20061 ) 



6. Summary and Outlook 

The canonical abundance of ^®A1 in the Solar 
System cannot be explained in terms of a late in- 
jection of debris from a nearby SN into the gaseous 
protoplanetary disk because (i) the dispersal of the 
natal cluster and the finite time window for injec- 
tion make it a priori an unlikely event (<2%), (ii) 
^^Al was already present in CAIs, which formed 
within ^ 10^ yr of the initial collapse of the pro- 
tosolar nebula and the formation of the protoplan- 
etary disk, and (Hi) the oxygen isotope systemat- 
ics of primitive Solar System materials show no 
sign of a late introduction of SN ejecta. The ap- 
parently uniform distribution of Al in meteorites 
and samples of the Earth, Moon, and Mars sug- 
gests homogenization during the collapse of the 
protosolar cloud. 
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We showed that the canonical ^^A1/^''A1 ratio 
can be explained if the Solar System formed from 
a molecular cloud contaminated by Wolf-Rayet 
winds from massive stars that formed 4-5 Myr ear- 
lier. A SN contribution is not required to explain 
the abundance of ^^Al, although it is not nec- 
essarily excluded. The a priori probability that 
such a level of contamination occured depends on 
the poorly-understood star formation histories in 
GMCs; we estimate that it is ~6%. However, our 
model predicts that the canonical value is close 
to the most likely non-zero value in the Galactic 
distribution. 

The initial "'^Ca/'*°Ca ratio can also be ex- 
plained by Wolf-Rayet wind contamination if 
~0.5 Myr elapsed between its introduction by 
winds and the formation of CAIs. If this scenario 
is also to explain primordial ^'^Fe in the Solar 
System, the cloud must have been contaminated 
with SN ejecta as well. If SN ejecta is included, 
our model over-predicts the abundance of ^°Fe by 
an order of magnitude. This discrepancy could be 
rectified by some combination of the following: (a) 
most ^°Fe falls back onto SN remnants rather than 



be ejected (jMever 2005t iTakigawa et al. 20081) : 

(b) most dust grains in SN ejecta are retained and 
destroyed in the shocked ejecta and never enter 
the molecular cloud; (c) the protosolar cloud was 
already collapsing and presented a smaller cross- 
section when the SN ejecta arrived; and (d) the 
^"Fe is a relict of an even earlier episode of mas- 
sive star formation and contamination for which 
all the other SLRs have decayed away. Absence 
of a significant excess of ^'^Mn seems to require 
(a), but not to the exclusion of the other expla- 
nations. Our model does not explain the inferred 
abundance of '^^Cl, and another mechanism such 
as irradiation much be invoked. 

A key uncertainty in our model is the efhciency 
with which ^^Al is introduced into the host molec- 
ular cloud and the degree to which it becomes uni- 
formly mixed. We propose that the carrier of Al 
was dust grains and that these dynamically decou- 
pled from the wind and embedded themselves (in- 
tact) into the cloud, but this hypothesis needs fur- 
ther investigation. Furthermore, our model does 
not account for the inhomogeneities in SN ejecta 
and WR winds that could produce spatial vari- 
ation in the contamination of a molecular cloud. 
There are also uncertainties in calculations of the 



evolution and nucleosynthesis of massive stars that 
could quantitatively alter our results. Production 
of ^^Al by neon burning during the Type Ib/c 
SN that follows the Wolf-R ayet phase is sensi - 
tive to the progenitor mass teig don et al. 20o3 ) 
and for a 60 M© progen itor, could be as large a s 
the yield from the wind (jWoosley fc Heger 20071 ). 
New models that include stellar rotation predict 
higher yields of ^^Al, an earlier appearance of Al 
in the WR wind (as early as 1 Myr) , and a smaller 
minimum initial m ass for entry into the WR phase 
([Palacios et al. 20 051. Larger ^^Al yields would 
relieve the requirement for high delivery efficiency 
to the molecular cloud and may resolve the dis- 
crepancy between the predicted amount of con- 
comitant ^"Fe from SN ejecta and the inferred ini- 
tial abundance of the radionuclide in the Solar Sys- 
tem. Future tests of this model could compare pre- 
dicted WR stellar contamination with short-lived 
isotopes (e.g., ^"""Pd) whose abundances seem con- 
sistent with models of the interstellar medium (G. 
Huss, pers. comm.), as well as the Solar System's 
oxygen isotopic composition. 
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Table 1 

Statistics of ^^Al abundance in different scenarios 



Scenario 


% 


26A1/27A1 




> > 5 X 10-5 


95% 


Disk injection: 






T = 0.5 Myr 


16 1.2 


5 X 10-6 


T = 1 Myr 


20 1.9 


1.3 X 10-5 


Disk injection with Bondi-Hoyle accretion 


T = 0.5 Myr 


16 1.2 


8 X lO-'^ 


T = 1 Myr 


21 1.9 


1.6 X 10-5 


Molecular gas contamination (T=0): 




2.7 Myr exp. SF 


16 6.2 


9 X 10-5 


1.7 Myr cxp. SF 


8 4.3 


6 X 10-5 


3.7 Myr exp. SF 


21 6.5 


9 X 10-5 



Note. — Only simulations which produced °A1- 
contaminated systems are reported. 
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0.15 





(c) . \ (d) . ' ; 

Fig. 1. — Cartoon of the molecular cloud contam- 
ination scenario to explain the presence of ^^Al 
in the early Solar System: (a) A young cluster 
of stars forms an HII region in a giant molecular 
cloud; (b) massive stars eject ^^Al into the cloud; 

(c) a clump collapses from the contaminated cloud; 

(d) a second generation of stars (including the 
Sun) forms from the clump. 



Fig. 2. — The distribution of non-zero values of the 
^^Al/^'^Al ratio predicted by a Monte Carlo model 
of its formation in massive stars and incorporation 
injection into the protoplanetary disk (see text for 
details). The units of the ordinate are fractional 
number of Monte Carlo systems per unit common 
logarithm. The vertical bar demarks the canonical 
Solar System value of 5 x 10~^. The curves do 
not integrate to unity because a large majority of 
systems are not contaminated with ^^Al (Tabled]). 
The solid curves are for the disk injection scenario 
(see text), and the dashed curves include Bondi- 
Hoyle accretion of nebular gas. Two values for the 
elapsed time between the formation of the Sun and 
^^Al-containing CAIs are considered. 
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Fig. 3. — The distribution of non-zero values of the 
^^Al/^'^Al ratio predicted by an alternative model 
in which ^^Al was produced in an earlier gener- 
ation of massive stars and introduced by Wolf- 
Rayet winds into the molecular cloud that formed 
the Sun (see text for details). The units of the 
ordinate are fractional number of Monte Carlo 
systems per unit common logarithm. The ver- 
tical bar is the canonical Solar System value of 
5 X 10~^. The integral of the curves are not unity 
because the majority of systems are not contam- 
inated with ^"^Al (Table [1]). The broken curves 
are for a monotonic elapsed time between the for- 
mation of the earlier generation of massive, ^^Al- 
producing stars and the Sun. The solid curves 
(shown one-tenth scale) are for an exponentially- 
decaying rate of star formation with a decay time 
of 2.7 Myr (heavy curve) or 1.7 and 3.7 Myr (light 
curves). The area under the solid curves changes 
with assume star formation history, but the loca- 
tion of the peak does not. 




'AI/"AI 



Fig. 4. — Calculated abundances of ^®A1 and 
^"Fe relative to stable comparison isotopes in 
star-forming regions contaminated by Wolf-Rayet 
winds and SN ejecta (see text for details). Each 
point represents a Monte Carlo calculation of the 
composition of the gas in a well-mixed molecular 
cloud 3-8 Myr after massive star fromation. The 
large black point is the inferred composition of the 
Solar System. The line represents the Galactic 
average abundance ratio from 7-ray observations 
(|Wang et al. 20071) . 
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